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ABSTRACT 

In this paper, we focus on the limb-brightened, prolate spheroidal cavity of the radio galaxy Cygnus A, 
as revealed by the Chandra X-ray Observatory. We use the shock heated, thermal intracluster medium 
around the expanding cavity to infer the properties of the radio synchrotron-emitting gas inside the 
cavity. The gas along the N and S edges of the cavity is found to have an average temperature of 6.0 
keV, which is hotter than the temperature (4.6 keV) of the adjacent intracluster gas. It is proposed 
that this hotter gas is intracluster gas shocked by the expanding cavity. The shock is thus inferred 
to be weak (Mach number 1.3, a value also inferred from the density jump at the cavity edge) and 
its velocity 1,460 km s" 1 . The total kinetic power of the expansion is found to be 1.2 x 10 erg 
s _1 , which is larger than both the total radio power and the power emitted by the entire intracluster 
medium in the 2-10 keV band. It appears that most of the power of the jets in Cygnus A is 
currently going into heating the intracluster medium. From the derived pressure inside the cavity, 
there is no conclusive evidence for a component contributing pressure additional to the magnetic fields 
and relativistic particles responsible for the synchrotron radio emission. Further, the ratio of energy 
densities in positive to negative cosmic rays in Cygnus A is between 1 and 100 (the value in our 
Galaxy). 

Subject headings: galaxies: active - galaxies: clusters - galaxies: individual (Cygnus A) - galaxies: jets 
- galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

It is now generally accepted that the lobes of powerful 
radio galaxies are fuelled by two oppositely directed rel- 
ativistic beams emitted from the nucleus. By comparing 
the power supplied by the beams with the power radi- 
ated by the radio components, Scheuer (1974) showed 
that only a small fraction of the power supplied can be 
radiated away from the tip of the beam, and that most 
of the energy must be deposited in a cavity surrounding 
the beam. The dynamics of this cavity (or "cocoon" as it 
is sometimes called) and its interaction with the external 
medium have been extensively discussed (e.g. Scheuer 
1974; Begelman & Cioffi 1989; Kaiser & Alexander 1997, 
1999; Heinz, Reynolds & Begelman 1998; Alexander 
2000; Reynolds, Heinz & Begelman 2001). 

In order to derive basic properties of the cavity, such 
as its pressure and expansion velocity, previous workers 
have had to use the synchrotron radio emission, often 
assuming equipartition of energy between magnetic fields 
and relativistic particles. There is no guarantee that this 
is the case and there may be other sources of pressure in 
the cavity, such as low energy cosmic rays or diffuse hot 
gas. A more direct probe of the properties of the cavity 
could be achieved if thermal emission from the shock 
driven by the overpressured cavity into the surrounding 
gas could be observed. 

In this paper, we deduce the presence and properties of 
the cavity and its external (bow) shock from Chandra X- 
ray observations of the thermal emission associated with 
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the powerful radio galaxy Cygnus A. Deficits of X-ray 
emission coincident with the inner radio lobes of Cyg A 
were first found by Carilli, Perley & Harris (1994) from 
ROSAT observations. However, ROSAT did not have 
the spatial or spectral resolution to derive the detailed 
morphology or gas temperatures in the vicinity of the 
cavity. The presence of a bow shock ahead of radio hot 
spot B of Cyg A was inferred by Carilli, Perley & Drchcr 
(1988) from measurements of Faraday rotation. 

In the present paper, Section 2 reviews the Chandra 
observations, while Section 3 derives various properties 
of the expanding cavity and the shock driven into the 
intracluster medium. Section 4 summarizes our conclu- 
sions. 

2. CHANDRA OBSERVATIONS 

The Chandra observations of Cyg A have been de- 
scribed in our earlier papers (Wilson et al. 2000; Young 
et al. 2002; Smith et al. 2002, hereafter Paper III). For 
consistency with most recent discussions of the topic of 
this paper, we use Ho = 75 km s _1 Mpc -1 and qo = 
(unlike Papers I - III which assumed Ho = 50 km s _1 
Mpc -1 and qo = 0). The luminosity distance to Cyg A 
is then dL = 231 Mpc and the angular size distance dA 
= 207 Mpc, so 1" corresponds to 1.0 kpc. Preliminary 
remarks on Chandra observations of the cavity have been 
presented in conference proceedings (Wilson et al. 2002, 
2003). 

For illustration purposes, Fig. 1 shows the Chandra 
image of Cyg A. The elliptical shaped morphology is in- 
terpreted in terms of hot gas swept up by the cavity. Fig. 
2 shows 6 cm VLA radio contours superposed on the X- 
ray grey scale. It can be seen that the radio lobes are as- 
sociated with relatively faint X-ray emission and project 
inside the limb-brightened edge of the cavity, both of 
which are in accord with expectation. 
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keV given in Table 4 of paper III. 

Fig. 3 shows the cavity with areas emitting bright 
X-ray emission outlined. The updated temperatures 
of MEKAL models are indicated. These regions are 
so bright that projection of foreground and background 
cluster emission onto them is negligible. In Fig. 3, re- 
gions with temperature > 5 keV are outlined, and the 
temperature labelled, in blue, while red is used for re- 
gions with temperatures < 5 keV. It is immediately note- 
able that the four regions projecting along the N and S 
edges of the cavity are hot (5.2 - 6.8 keV), while those 
in the "belts" crossing the cavity and the nucleus of the 
galaxy from NE to SW are cool (3.8 - 4.4 keV). The 
"belts" (which are not discussed in this paper) appear 
to be a disklike structure of cooling gas, seen obliquely 
and being accreted into the galaxy. This cooling gas may 
form the outer part of the accretion disk which fuels the 
nuclear activity. 
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Fig. 2. — Contours of a VLA 6 cm radio map (Perley et al. 1984) 
superposed on a grey scale image of the X-ray emission. 



3.1. Mach Number and Shock Velocity of Expanding 

Cavity 

We interpret the hot regions elongated along the N and 
S edges of the cavity as intracluster gas shocked by the 
expanding cavity. The fact that these hot regions are 
visible at all implies that the stand-off distance of the 
shock is comparable to or larger than Chandra's resolu- 
tion (about 1 kpc). The deprojected temperature of the 
innermost shell of the cluster, T\ = 4.60 keV, is taken as 
the pre-shock temperature. The average temperature - 
T-x = 6.01 keV - of the four "blue" regions on the N and 
S edges of the cavity is assumed to be representative of 
the postshock temperature. 

For a nonradiating shock, equations (10-19) and (10- 
20) of Spitzer (1978) are 

P2 _ 2jM_ 2 7 - 1 
Pi 7 + 1 7+1 

and 



3. EXPANSION OF THE CAVITY 

In Paper III, we deprojected the X-ray spectra of the 
intracluster gas around Cygnus A taken from 12 elliptical 
and circular annuli outside the cavity in order to derive 
the run of temperature, metallicity, density and pressure 
as a function of radius. The goal was to derive the ambi- 
ent values of the parameters presumed unaffected by the 
energy input by Cygnus A itself. The temperature of the 
X-ray emitting gas was found to drop from ~ 8 keV more 
than ~ 80 kpc from the center to ~ 5 keV some 50 kpc 
from the center, with the coolest gas immediately out- 
side the cavity. Since this work was published, improved 
response functions ("fefs") have become available. We 
have, therefore, deprojected each cluster shell and mod- 
elled each one with a single temperature MEKAL plus 
foreground emission, absorbed by a column of cold gas, 
precisely as done in Paper III but using the new fefs. 
Such a model provides an acceptable fit to the spectra of 
all 12 shells, and gives small changes in the run of param- 
eters with radius compared with the results presented in 
Paper III. In particular, the (deprojected) temperature 
in the innermost shell, just outside the cavity, is now T\ 
= 4.60 ± 0.5 keV, consistent with the value of 4.91 ± 0.6 



m = Ei = i-i + f ( 2 ) 

where pi {p 2 ), «i (1*2) and p\ (p 2 ) represent the pre-shock 
(post-shock) pressure, velocity and density, 7 is the ratio 
of the specific heats and M the Mach number. We have 
assumed that the magnetic field exterior to the cavity 
is sufficiently weak that it does not affect the dynamics. 
From these equations, the ratio of the post-shock (T 2 ) to 
pre-shock (2\) temperatures is 

T 2 [2 7 M 2 -( 7 -l)][(7-l) + (2/M 2 )] _ 
T X ~ (7 + 1) 2 ■ [ ' 

Equation (3) can be inverted to a quadratic equation in 
M 2 . Putting 7 = 5/3, T 2 = 6.01 keV, and Ti = 4.60 keV, 
we thus find the Mach number of the shock driven by the 
expanding cavity to be M = 1.32. We note in passing 
that Alexander & Pooley (1996) also concluded that the 
hot spot advance speed is only mildly supersonic (M ~ 
2) on completely different grounds. 

The thermal bremsstrahlung emissivity depends only 
weakly on temperature (e„ cx T -1 / 2 ). We can thus use 
the jump in emissivity at the limb-brightened edge of the 
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Fig. 3. — An unsmoothed image of the central region of Cyg A in the 0.75 - 8 keV band. The areas indicated by solid lines mark regions 
from which spectra of the X-ray emission were extracted and modelled with MEKALs; the numbers indicate the resulting gas temperature 
in keV. Temperatures of regions marked in blue are hotter than the innermost cluster shell (4.60 keV), while those marked in red are cooler. 
It is noteable that the hotter regions are associated with the limb-brightened edges of the cavity, while the cooler regions are in the "belts" 
which traverse or encircle the cavity along its minor axis, and may represent the outer regions of the accretion disk. 



cavity (a factor of ~ 3.1 in p. a. 30° - 40° from the nu- 
cleus) to make a semi- independent estimate of M. We 
find a density ratio p\f pi — 0.67, and a Mach number 
M ~ 1.3 (equation 2). However, within the uncertainties 
of both methods, it is possible that Cyg A is simply ex- 
panding at the sound speed of the external intracluster 
gas. We shall, nevertheless, continue to assume M = 1.3 
in the calculations that follow. Both methods have their 
uncertainties from measurements, deprojections and the 
simple analytical model, but the excellent agreement be- 
tween the two is encouraging. 

The shock velocity may now be found from the relation 
(Spitzer 1978, equation 10-21) 



The total kinetic power of the expanding cavity is 



1 kT 1 M 2 



(4) 



where p is the mean mass per particle and k is the Boltz- 
mann constant. Assuming the above values for 7, T\ and 
M, and that p, = 1.0 X 10" 24 gm, we find ui = 1,460 
km s" 1 . All values derived so far are independent of 
the Hubble constant. Because we have used gas with a 
range of azimuthal angles around the cavity to derive this 
shock velocity (at least in the "temperature method"), 
we will take u\ = 1,460 km s _1 as an approximation to 
the angle-averaged velocity of the shock driven by the 
expanding cavity. 

3.2. Luminosities, Timescales and Pressures 



L. 



1 



-piu\A : 



(5) 



where A is the total area of the shock surface. We assume 
this surface has a prolate spheroidal shape, with major 
axis in the plane of the sky. From the X-ray image, we 
adopt a semi-major axis of a = 62" and a semi-minor 
axis of b = 28" and find A = 1.7 x 10 47 cm 2 . From 
the deprojection of the cluster emission, the density of 
the innermost cluster shell, which we take to be the pre- 
shock density, is p\ — 4.4 x 10" 26 gm cm" 3 . With u\ = 
1,460 km s _1 , we have L exp ~ 1.2 x 10 46 erg s _1 . For a 
black hole mass of (2.5 ± 0.7) x 10 9 M (Tadhunter et 
al. 2003), L exp /L Edd is 0.037. 

L exp exceeds the total radio emission of Cyg A (Lr(10 
MHz - 400 GHz) = 9.0 x 10 44 erg s" 1 ) and the total clus- 
ter gas X-ray emission in the 2 to 10 keV band {Lx(2 
- 10 keV) = 1.2 x 10 45 erg s _1 . These values are from 
Carilli & Barthel (1996) adjusted to our assumed lumi- 
nosity distance. The dominant current power source in 
Cyg A thus appears to be the expansion of the jet-driven 
cavity which heats the surrounding intracluster gas. L exp 
represents a lower limit to the current power of the two 
jets. 

We assume the cavity maintains its axial ratio (a/6 ~ 
2.2) as it expands. The shock velocity along the minor 
axis would then be approximately Ub — (b/a 
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Ui 



4 



The Cavity of Cygnus A 



992 km s _1 and along the major axis u a ~ [a/b) 1 ' 2 u\ 
= 2,150 km s -1 . Assuming these velocities represent 
the velocities of expansion of the cavity and have been 
constant over the lifetime of the radio source, the age of 
Cyg A is t ago ~ 3.0 x 10 7 yrs. This value is within a 
factor of two of other estimates of the age based on the 
radio properties (e.g. Begelman & Cioffi 1989; Carilli ct 
al. 1991; Kaiser & Alexander 1999). We emphasise that 
any net infall ("cooling flow" or gravitational infall) of 
the pre-shock gas has been assumed to be zero here, as 
we have no way of determining it. 

The cooling timescale of the hot, postshock gas, with 
p 2 ~ 6.5 x 10~ 26 gm cm- 3 and T 2 ~ 7.0 x 10 7 K, is t coo i 
~ 2.5 x 10 8 yrs. This is a factor of 8 times larger than 
the age of Cyg A. Thus, the expanding cavity is currently 
heating up the surrounding gas; cooling of the gas is a 
slower process. This finding that t coo i > t ago supports 
our assumption of a nonradiating shock. 

The lateral expansion of the cavity is driven by its 
internal pressure (Scheuer 1974): 

Piul = Pc (6) 

With the values given above for p\ and Ub, we have p c ~ 
4.4 x 10~ 10 dyne cm~ 2 . If this pressure is produced by 
(tangled) magnetic fields and isotropic relativistic parti- 
cles (so p c = u c /3) in cquipartition, the magnetic field in 
the cavity is B c ~ 1.2 x 10~ 4 gauss. This value is close 
to the equipartition magnetic field strength in the lobes 
calculated from their radio emission (6 x 10~ J gauss for 
our distance, Winter et al. 1980; (4.5 - 6.5) x 10" 5 
gauss, Carilli et al. 1991; (5-7) x 10~ 5 gauss, Alexan- 
der & Pooley 1996). The discrepancy between our value 
of B c and that obtained from the radio observations may 
not be as significant as it seems given the uncertainties 
in both estimates. 

We note, in addition, that all the equipartition field 
strengths in the lobes quoted above assumed equal en- 
ergy densities in positive particles and electrons (C. Car- 
illi, P. Alexander, private communications). If it is as- 
sumed, instead, that the positive particles have 100 times 
the energy density of the electrons (as is the case in our 
Galaxy), then the equipartition field strengths would be 
3.7 times larger and the pressure in cosmic rays and mag- 
netic fields 14 times larger. Such a pressure is larger than 
deduced above from equation (6), so we do not favor the 
notion that the energy density in positive particles is 100 
times that in electrons, but we cannot rule it out. 

We have found the shock velocity along the major axis 
of the cavity to be u a ~ 2,150 km s _1 . This agrees well 
with the hot spot advance speed of 1,500 km s _1 found 
by Alexander & Pooley (1996), assuming ram pressure 
confinement and that the whole front surface of the radio 
source moves forward with the same speed. 



4. CONCLUSIONS 

In this paper, we have explored the dynamics of the 
cavity (or "cocoon") associated with Cygnus A. The pri- 
mary source of information was Chandra X-ray imaging 
and spectroscopy, from which densities and temperatures 
of the thermal gas in the intracluster gas around the cav- 
ity were obtained. 

Our primary assumption is that the hotter (6.0 keV) 
gas in the limb-brightened regions at the N and S edges 
of the cavity has been shocked by the expanding cavity. 
From this assumption, we infer the Mach number (1.3) 
and velocity (1,460 km s" 1 ) of the shock. The same num- 
bers are found from the density jump. If, instead, some 
of this gas has not been shocked, we will have underesti- 
mated the postshock temperature and thus also the Mach 
number, expansion velocity and power of the expanding 
cavity. The kinetic power of the expanding cavity (~ 1.2 
x 10 46 erg s _1 ) exceeds the total radio emission in the 10 
MHz to 400 GHz band by a factor of ~ 13 and the total 
X-ray emission from the cluster in the 2-10 keV band 
by a factor of ~ 10. The cooling time of the shocked gas 
(2.5 x 10 s yrs) is longer than the age of the radio source 
by a factor of ~ 8. It thus appears that most of the 
power currently emitted in the relativistic beams which 
fuel the radio source goes into heating the intracluster 
medium. We find that the pressure inside the cavity (4.4 
x 10~ 10 dynes cm~ 2 ) exceeds the minimum pressure in 
magnetic fields and cosmic rays in the radio lobes. This 
difference may not be significant, so we find no conclu- 
sive evidence for any other entity contributing pressure 
inside the cavity. A further consequence is that the ratio 
of the energy density in positive relativistic particles to 
that in relativistic electrons in Cyg A is less than the 
value of 100 found in our Galaxy. 

If the cavity shock is weak (M = 1.3) in a powerful, 
young radio galaxy like Cyg A, then it seems that, in gen- 
eral, heating of the intracluster medium by radio galax- 
ies must be through the media of weak shocks and sound 
waves. Strong shocks are implausible (cf. Heinz 2003). 

The accuracy of our results is limited by the simpli- 
fied analytical model used. Further work, combining the 
results of the X-ray observations with modern computa- 
tional simulations of Cyg A, would be well worthwhile. 
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